Cross-Correlation Studies with CMB Polarization Maps 
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The free-electron population during the reionized epoch rescatters CMB temperature quadrupole 
and generates a now well-known polarization signal at large angular scales. While this contribution 
has been detected in the temperature-polarization cross power spectrum measured with WMAP 
data, due to the large cosmic variance associated with anisotropy measurements at tens of degree 
angular scales only limited information related to reionization, such as the optical depth to electron 
scattering, can be extracted. The inhomogeneities in the free-electron population lead to an ad- 
ditional secondary polarization anisotropy contribution at arcminute scales. While the fluctuation 
amplitude, relative to dominant primordial fluctuations, is small, we suggest that a cross-correlation 
between arcminute scale CMB polarization data and a tracer fleld of the high redshift universe, 
such as through fluctuations captured by the 21 cm neutral Hydrogen background or those in the 
infrared background related to flrst proto-galaxies, may allow one to study additional details related 
to reionization. For this purpose, we discuss an optimized higher order correlation measurement, in 
the form of a three-point function, including information from large angular scale CMB temperature 
anisotropics in addition to arcminute scale polarization signal related to inhomogeneous reionization. 
We suggest that the proposed bispectrum can be measured with a substantial signal-to-noise ratio 
and does not require all-sky maps of CMB polarization or that of the tracer field. A measurement 
such as the one proposed may allow one to establish the epoch when CMB polarization related to 
reionization is generated and to address if the universe was reionized once or twice. 



I. INTRODUCTION 

The increase in sensitivity of upcoming cosmic mi- 
crowave background (CMB) polarization experiments, 
both from ground and space, raises the possibility for 
detailed studies related to reionization. The main rea- 
son for this is the existence of a large angular scale po- 
larization contribution due to rescattering of the tem- 
perature quadrupole by free-electrons in the reionized 
epoch [1] . This contribution peaks at angular scales cor- 
responding to the horizon at the rescattering surface. 
While such a signal has now been detected with the 
temperature-polarization cross-correlation power spec- 
trum measured with first-year data of the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) [2], information re- 
lated to reionization from it, unfortunately, is limited. 
While the amplitude of the signal depends on the total 
optical depth to electron scattering, with an estimated 
value of 0.17 ± 0.04 [3], the reionization history related 
to this optical depth is still unknown. While slight mod- 
ifications to the large angular scale polarization power 
spectra exist with complex reionization histories, such as 
a two stage reionization process advocated by Ref. [4], 
due to large cosmic variance associated with anisotropy 
measurements at few tens degree scales, one cannot fully 
reconstruct the reionization history as a function of red- 
shift [5]. 

Given the importance of understanding reionization for 
both cosmological and astrophysical reasons, and the in- 
crease in sensitivity and angular resolution of upcoming 
CMB polarization data, it is useful to consider alterna- 
tive possibilities for additional studies beyond the large 
angular scale signal. Under standard expectations for 
reionization, mainly due to UV light emitted by first lu- 



minous objects, the process of reionization is expected 
to be patchy and inhomogeneous [6]. This leads to fluc- 
tuations in the electron scattering optical depth and a 
modulation to the polarization contribution such that 
new anisotropy fluctuations, at arcminute scales corre- 
sponding to inhomogeneities in the visibility function, is 
generated [7,8]. Even if the reionization were to be more 
uniform, as expected from alternative models involving 
X-ray backgrounds [9] and reionization via particle de- 
cays [10], density fluctuations in the electron field can 
modulate the polarization contribution and generate sec- 
ondary anisotropics [11]. The new secondary polariza- 
tion fluctuations, whether due to patchiness or density 
inhomogeneities, however, is hardly detectable in the po- 
larization power spectrum given the dominant primor- 
dial polarization contribution, with a peak at arcminute 
scales, related to the velocity field at the last scatter- 
ing surface. Thus, information related to reionization is 
limited from measurements that involve simply the two 
point correlation function of polarization or the cross- 
correlation between polarization and the temperature. 
To extract additional information one must move to a 
higher order and consider measurements related to, for 
example, the three-point correlation function or the bis- 
pectrum [12]. 

In this paper, we consider a potentially interesting 
study that has the capability to allow detailed mea- 
surements related to the reionization history. The pro- 
posed measurement involves the cross-correlation be- 
tween CMB polarization maps and images of the high 
redshift universe around the era of reionization or its end. 
The existence of such a correlation arises from the fact 
that the polarization contribution related to reionization 
is generated at the same epoch as that related to the 
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tracer field selected from high redshifts. Since the reion- 
ization contribution to polarization is generated by the 
scattering of the temperature quadrupole, direct cross- 
correlation between polarization maps and a tracer field 
is not useful. However, since the rescattering quadrupole 
correlates with CMB temperatiirc maps, one must con- 
sider higher order statistics that also use information 
from CMB temperature data. We consider a possibil- 
ity in the form of a three point correlation function, or 
a bispectrum in Fourier space, involving CMB polariza- 
tion maps at arcminute-scale resolution, CMB tempera- 
ture maps with large angular scale anisotropy informa- 
tion, and maps of the high redshift universe that trace 
fluctuations during the reionization and prior to reion- 
ization. Examples of useful fields include maps of the in- 
frared background anisotropics related to the redshifted 
UV emission by first proto galaxies [13] and fluctuations 
in the 21 cm background due to the neutral Hydrogen 
content [14]. Additional possibilities also include maps of 
the 2^3 universe since one can then establish the low- 
redshift part of the electron scattering visibility function. 
Cross-correlation studies with the 21 cm background may 
be the most interesting possibility given that one is al- 
lowed to a priori select redshift ranges from which fluc- 
tuations arise to the redshifted 21 cm emission, based on 
frequency information of the observations. 

The paper is organized as follows. In § H, we derive 
the existence of a higher order correlation between CMB 
temperature, polarization and a tracer field of the large 
scale structure that overlaps with fluctuations in the visi- 
bility hmction. In § III, we discuss our results and suggest 
that there is adequate signal-to-noise to perform such a 
measurements and consider several applications. We con- 
clude with a summary in § IV. 

II. CALCULATION METHOD 

The rescattering of CMB photons by free-electrons at 
redshifts significantly below decoupling, at z ^ 1100, 
leads to a polarized intensity, which in terms of the linear 
Stokes parameters, is 

±X(n) = 

I 2 

{^±^u)(r.) = ^ fdrg{r) ^ g(™) ±2!^ (A) , W 

m=-2 

where the quadrupole anisotropy of the radiation field is 

dO-|=e(x,n). (2) 

In Eq. 1, r{z) is the conformal distance (or lookback time) 
from the observer at redshift z = 0, given by 



where the expansion rate for adiabatic CDM cosmological 
models with a cosmological constant is 

= Hi [17^(1 + zf + nK{l + zf + 17a] , (4) 

where Hq can be written as the inverse Hubble distance 
today = 2997.9ft,"iMpc. We follow the conventions 
that in units of the critical density STJq/SttG, the con- 
tribution of each component is denoted 17 j, i = c for the 
CDM, h for the baryons, A for the cosmological constant. 
We also define the auxiliary quantities f7m = fJc + and 
0,K = 1 — 17i, which represent the matter density and 
the contribution of spatial curvature to the expansion 
rate respectively. Although we maintain generality in all 
derivations, we illustrate our results with the currently 
favored ACDM cosmological model. The parameters for 
this model are f2c = 0.30, f2(, = 0.05, Oa = 0.65 and 
h = 0.65. 

The visibility function, or the probability of scattering 
within dr of r, is 

g = Te-^ = X^{z)HQTH{\ + zfe-^. (5) 

Here r(r) = drf is the optical depth out to r, Xf.{z) 
is the ionization fraction, as a function of redshift, and 

TK = 0.0691(1 - Yp)Q.bh , (6) 

is the optical depth to Thomson scattering to the Hubble 
distance today, assuming full hydrogen ionization with 
primordial helium fraction of Yp{= 0.24). In typical cal- 
culations of CMB anisotropics, the general assumption 
is that the universe reionized smoothly and promptly 
such that Xe{z < Zri) = 1 and Xe{z > Zri) = 0. Such 
an assumption, however, is in conflict with observations, 
so far, involving the optical depth to electron scatter- 
ing measured by WMAP and the presence of a 1% neu- 
tral fraction based on Lyman-a optical depths related 
to Gunn-Peterson troughs of the z ^ 6 quasars in the 
Sloan Digital Sky Survey [15]. To explain the WMAP 
optical depth and low redshift data simultaneously re- 
quire complex reionization histories where Xe{z) is not 
abrupt, but where reionization takes a long time [16]. 
Since first luminous objects are usually considered as a 
source of reionization, the process is inhomogeneous and 
patchy. 

Here, we make use of two descriptions of Xe{z) (Fig. 1; 
left panel). In the first case (model A), we use a calcula- 
tion based on the Press-Schechter [17] mass function and 
related to the reionization by UV light from the first-star 
formation following Rcf. [18]. Here, Xe{z) varies from 
a value less than 10~^ at a redshift of 30 to a value of 
unity when the universe is fully reionized, at a redshift 
of ~ 5. The reionization history is rather broad and the 
universe does not become fully reionized till late times, 
though the reionization process began at a much higher 
epoch. The total optical depth related to this ionization 
history is ~ 0.17, consistent with WMAP measurements 
[3]. As an alternative to such a history, and to consider 
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FIG. 1. Le]t -panel: The reionization history of the universe. Here, we consider two possibihties for the reionization process 
involving a smooth transition over a wide range of redshift (model A) and a model which involve a sudden transition to a 
fully reionized model at a reionization redshift ~ 15. Right panel: The visibility function related to electron scattering in the 
rcionized epoch, g(z) (see, equation 5) for both models A and B shown above. For comparison, we also show several normalized 
redshift distributions of tracer fields that wc consider here. These are, (a) a high redshift broad distribution, such as the one 
related to first proto-galaxy contributions to the IR background, (b) a high redshift, but narrow, distribution expected from 
the 21cm background fluctuation studies in the future, and (c) a low redshift distribution, such as the distribution of « ~ 3 
galaxy population. Note that 21 cm fluctuation studies allow one to choose narrow bands, such as (b), in redshift over a wide 
range both prior to and during the reionization process and may provide the optimal tracer fleld for the proposed study here. 



the possibility how different reionization models can be 
studied with the proposed cross-correlation analysis, we 
also consider a reionization history that is instantaneous 
at the redshift of reionization of ~ 15 (model B). This 
alternative model could arise from reionization descrip- 
tions that involve the presence of a X-ray background or 
where the reionization process is associated with decay- 
ing particles, among others. While the transition to a 
fully reionized universe, from a neutral one, is sudden, a 
history such as the one shown produces a similar optical 
depth to electron scattering as in model A. 

In Fig. 1 right panel, we show the visibility function 
g{z) related to these two reionization histories, in model 
A, scattering is distributed widely given the long reion- 
ization process, while in the case of model B, scattering 
is mostly concentrated during the sharp transition to a 
reionized universe. As shown in Fig. 1, one does not ex- 
pect a probe of the high-z universe, say at a redshift ~ 20, 
to correlate strongly with polarization if the reionization 
process is more like model B, while a strong correlation 
is expected if the hypothesis related to a lengthy reion- 
ization process does in fact happen. 

Now we will discuss how these reionization histories 
determine polarization signals in CMB both at large and 
small angular scales. First, wc note that the polarization 
is a spin-2 field and can be decomposed using the spin- 
spherical harmonics [19] such that 

±X(n) = J2 ±Xi^ ±2Yr{i^) • (7) 

Im 



The multipole moments of polarization field are generally 
separated to ones with gradient {E) and curl [B) parity 
such that [20] 

±Xlm = Elm ± iBlm ■ (8) 

Instead of Stokcs-Q and -U, one can redefine two polar- 
ization related fields -E'(n), a scalar, and B{n), a pseudo- 
scalar, such that 

Im 

B{n) = J2BlmYrin). (9) 

I'm 

The cross correlation between reionization generated 
polarization and a tracer field of the large scale structure 
is zero. This can be understood based on the fact that 

(±X(n)5(n')) oc ^^(^^^(fi')) = q , (10) 

since the temperature quadrupole that rescatters at low 
redshifts is not correlated with the local density field. 
There is, however, a non-zero cross-correlation between 
polarization and the CMB temperature anisotropy due 
to the fact that the temperature quadrupole related to 
scattering is present in the CMB temperature anisotropy 
map: 

(±X(n)e(n')) oc ^(g(™)e(n')) ^ . (11) 

m 
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Since this is part of the standard calculation related to 
the Cf^ angular power spectrum, we do not further con- 
sider it. 

To the next highest order, contributions to cross- 
correlation arise from the fact that fluctuations in the 
electron field is correlated with the large scale structure. 
This second-order term arises from the fact that 



5f(n,r) = 5f(n,r) 



1 + 



5g{n,r) 



5(n,r) 



(12) 



These fluctuations lead to a second order polarization 
signal that peaks at arcminute angular scales and are 
discussed in Rcf. [7,11] 

Substituting Eq. 12 in Eq. 1 and using Eq. 7, we can 
write the multipole moments the polarization field, in the 
presence of fluctuations in the visibility function, as 



xy™Mki)CMk2)±2rr-^(ki) 

X J dn ±2yr(ri)r™nft)>^r(n)±2>r^ , (13) 

where, following Ref. [7], we have projected the 
fluctuations-modulated temperature quadrupole 



Qf'(k) = J ^Q('")(ki)<5,(|k-ki|), 



(14) 



to a basis where the ^;-axis is parallel to the k vector such 
that 

Q^Hk) = y^y^g(0)(kr)<5,(k2)Fr(ki), (15) 

and have simplified using the Rayleigh expansion of a 
plane wave 

gik.flr ^ 4;r^i'i((A:r)y(""(k)r™(n) . (16) 

Im 

Since what correlates with the polarization field is the 
quadrupole generated by large angular scale temperature 
fluctuations, to simplify, we make the assumption that 
the relevant anisotropies are due to the Sachs- Wolfe (SW; 
[21]) effect 



e(x,n) = --$(x,n) 



(17) 



where ro = r{z = 1100). Note that there is a correction 
here related to the integrated Sachs- Wolfe (ISW) effect, 
but since it is present only at rcdshifts less than 1, when 
the dark energy component begins to dominate the en- 
ergy density of the universe, we ignore its contribution 



to the quadrupole. This is a valid assumption for models 
of reionization with an optical depth at the level of 0.1 
or more since most of the rescattering is then restricted 
to redshifts greater than 10 or more*. 

Our assumption related to large angular scale fluctu- 
ations allows us to write the quadrupole at late times, 
when projected to an observer at a distance r through 
free-streaming, as 

Q(°)(k,r) = -y5i$(k,ro)j2(kr,), (18) 

where = rg — r. 

Using Eq. 17, multipole moments of the temperature 
map at large angular scales are 

Qim = -^i' j ^$(k,ro)jKfcro)Fr(k). (19) 

Similarly, multipole moments of the tracer field are 

Sim = J dnYr*{n)S{h), (20) 

and assuming that the tracer field can be described with 
a source radial distribution of W^{r) with fluctuations 
denoted by 6s, we simplify to write 

Sirr,=4W [ I drW^{r)5s{\^,r)ji{kr)Yr{\^). 

J (27r) J 

(21) 

The cross-correlation with the large scale structure ex- 
ists in the form of a bispectrum between polarization- 
temperature-tracer fields. For this, we construct, for 
example, {Ei^m^Qi^m-iSi^mi) and after some straightfor- 
ward but tedious algebra, we write 

(-E'iimi 0i2m2 "^'sTOa ) — 

At, J dn J dr2lff{n,r2)JtM) 

X JdnYi;^^in)Y,^-^(h)Y,^^^{h), (22) 

where 

/f^(ri,r2) = ^J k^dkg{n)W^{r2)jk{kr2)ji,{kn)Pgs{k), 



Ifdk 
18^ 



P^t>{ki,ro)jiikro)j2{krs)ei{kr) . (23) 



Here, e/(x) = -ji{x)+j'^{x)+2ji{x)/x'^+Aj'i{x)/x and we 
have defined three-dimensional power spectra related to 



*When calculating the signal-to-noise, however, we include 

all contributions to anisotropics, such that the temperature 
anisotropy angular power spectrum, Cp'^^, is not just due to 
the SW effect. 
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the potential-field at the last scattering surface and the 
cross-power spectrum between fluctuations in the scat- 
tering visibility function and the tracer field: 



($(k)$(k')) = (27r)3,5°(k + k')P$$(fc) 
(<5,(k)5s(k')) -(27r)35D(k + k')P,5(fc), 



(24) 



respectively, where is the Dirac delta function. 

In these calculations, we will described Pgs{k,z) fol- 
lowing the halo model [22]. At large angular scales, 
relevant to most of the estimates in here, Pgs{k,z) = 
bg{z)bs{z)G'^{z)P{k), where P{k) is the power spectrum 
of density fluctuations in linear perturbation theory, and 
bg{z) and bs{z) are large physical scale biases of these 
two fields taken to be independent of the scale (see, be- 
low). We use the fitting formulae of Ref. [23] in eval- 
uating the transfer function for CDM models, while we 
adopt the COBE normalization to normalize the linear 
power spectrum [24]. This model has mass fluctuations 
on the 8h Mpc"""^ scale in accord with the abundance of 
galaxy clusters such that cg — 0.86. 

The bispectrum for the case with E-mode map is given 

by 



h h h 
mi 1712 ™3 



{-^hmi ^l2m2 ^tsms 



47r I j '^-'a ' 



(25) 



where. 



(26) 

One can simplify Eq. 26 with the use of the Limber 
approximation. Here, we use a version based on the com- 
pleteness relation of spherical Bessel functions [12] 

/ dkk'F{k)ji{kr)ji{kr') « 5d^'<5°(r - r')F{k)\^^^ , 

(27) 

where the assumption is that F{k) is a slowly- varying 
function and the angular diameter distance in above is 



dA = H^^n^^'"^ sinh(iJof^^ V) 



1/2 



(28) 



Note that as VIk ^ 0, dyi ^ r. 

Applying this approximation to the integral in If^ 
(Eq. 23) allows us to further simplify and write 



b^ - — 



drg^^W^{r)Pgs[k. 



dA 



iFAr), (29) 




FIG. 2. The integrals 7f (r), with thick lines (Eq. 30), and 
iP{r), with thin lines (Eq. 33), as a function of the multipole 
for several different redshift values. These mode coupling 
integrals oscillate between positive and negative values and 
we label positive/negative parts with (+)/(-). These func- 
tions captures the cross-correlation between density modu- 
lated temperature quadrupole that rescatter by an electron, 
at a radial distance of r, to form E- and B-modes of polar- 
ization and the temperature anisotropics seen by an observer 
at a redshift of zero or today. The multipole I here is that of 
temperature anisotropics as seen today. 



where 



iFir)^ / k^dkP^^{ki,ro)ji{kr„)j2{kr,)ei{kr). 



(30) 

Similar to the bispectrum involving the E-mode map, 
combined with temperature and a tracer field of the high- 
redshift universe, one can also construct a bispectrum in- 
volving the B-mode map instead of the E-mode map. As 
discussed in Ref. [7], rescattering under density fluctu- 
ations to the electron distribution also lead to a contri- 
bution to the B-mode map and can be understood based 
on the fact that the density-modulated quadrupole gener- 
ates m = ±1 components of the quadrupole which when 
scattered lead to B-modes^ . Following our earlier discus- 
sion, this bispectrum is written as 



E 



h h h 

mi 7712 ITT-S 



limi C;2'n2 ^hms 



(2;i + l)(2;2 + l)(2Z3 + l) 



47r 



h h h 




^Without density fluctuations, the temperature quadrupole 
related to scalar fluctuations induces only the m = compo- 
nent and only E-modes of the polarization are generated at 
large angular scales. 
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where, 
where 

7f(r)= J k''dkP^^{h,ro)ji{kro)j2ikr,)pi{kr) , (33) 

where, the function that replaces e; in to obtain ip 
ispi{x) = 2j'i{x)+ijl{x)/x. 

In Fig. 2, wc show the two quantities if'^ at several 
different values of z, as a function of I. While the function 
peaks at multipoles less than 4, at 2; = 1, the function 
broadens as one moves to rcdshifts greater than 10 such 
that it is basically a constant for / values up to ~ 20 or 
more. The quantities lf'^{r) only depends on the back- 
ground cosmology and the primordial power spectrum; 
with temperature anisotropy measurements and other 
cosmological studies, these parameters can assumed to 
be known. Thus, the only unknown qiiantity will be re- 
lated to reionization through the visibility fmiction, g{r). 
Thus, the cross-correlation between tracer fields allow 
some information related to g(r) be established. This in- 
cludes not only whether g{r) is inhomogeneous (through 
the detection of fluctuations), but also information on 
the form of g{r), as a function of r or binned r. 

III. RESULTS 
A. Cross-correlation 

It is now well known that the large scale tempera- 
ture fluctuations correlate with fluctuations in polariza- 
tion generated by rescattering due to free-electrons in 
the reionized epoch. Though there is no direct cross- 
correlation between polarization data and maps of the 
large scale structure, we suggest that the correlation 
between fluctuations in the large scale structure den- 
sity fields and that of the electron scattering visibil- 
ity function can be extracted by constructing a three- 
point correlation function, or a bispectrum in Fourier 
space, involving the large angular scale CMB temper- 
ature anisotropies, arcminutc scale polarization fluctua- 
tions, and maps of the high redshift universe. The re- 
quirements for the existence of such a correlation are (1) 
the presence of fluctuations in the visibility function, due 
to density inhomogeneities or patchiness of the reioniza- 
tion process, and (2) overlap between the redshift dis- 
tribution of scattering electrons, or the mean visibility 
function, and that of the tracer field. 

In Fig. 1, we illustrate the basis for the proposed study. 
Here, we show two visibility functions related electron 
scattering based on two different descriptions of reion- 
ization history consistent with current estimates of the 
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FIG. 3. The angular power spectrum of CMB polarizer 
tion in both the dominant gradient- or the E-mode and the 
B-modc. In the case of E-modcs, wc show the primordial con- 
tribution with the correction at large angular scales related 
to reionization scattering, and assuming a total optical depth 
to scattering of 0.17. In the case of B-modes, we neglect the 
primordial contribution related to gravity waves — both due 
to the unknown amplitude and the fact that this contribu- 
tion peaks at large angular scales while the proposed correla- 
tion is related to arcminute scales — and show the important 
arcminute-scale signal related to lensing conversion of E- to 
B-modes. The long-dashed and dot-dashed curves are the 
secondary polarization contributions, with power in E-modes 
equal to B-modes, related to inhomogeneities in the reioniza- 
tion involving ionized fraction and the electron density field, 
respectively. For reionization histories related to the forma- 
tion of first luminous objects, the reionization is significantly 
inhomogeneous and leads to a higher secondary contribution 
related to patchiness of the reionization when compared to 
fiuctuations in the electron density. The two contributions 
also arc generated at different epochs, with the patchy con- 
tribution resulting from redshifts where < Xe{z) < 1, while 
the density fluctuations generate contributions even when 
Xeiz) = 1. 

optical depth; In addition to these two, note that one 
can conceive a large number of reionization history mod- 
els that given the same optical depth [16]. In the first 
scenario involving reionization by the UV light from first 
stars (model A), the visibility function is rather broad, 
but generally peaks at redshifts where Xf.{z) ~ 0.5, while 
in the second scenario with a sudden transition to a reion- 
ized universe, scattering happens mostly during this tran- 
sition. 

To describe the three-dimensional cross power spec- 
trum between the tracer field and that of the visibility 
function, we make use of the halo model, but concentrate 
only on large angular scale clustering captured by the 2- 
halo term. In this limit, the relevant bias factors can be 
described as 
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FIG. 4. fcij^ig as a function of the mutlipole for the two 
reionization models separated to thick (model A) and thin 
(model B) lines. In the top panel, we show bi^^i^ as a func- 
tion of h for several different values of I3 while in the bottom 
panel we show bi^^i^ as a function of I3 for several values of 
h- The curves labeled (+)/(-) are positive/negative values of 
this quantity. Since bi^.i^ oc J drf{l-j,)Ii2 the top panel is sim- 
ply reflects Jj^, as shown in Fig. 2. Here, we assume a tracer 
field with a distribution given by the broad high-z distribu- 
tion labeled (a) in Fig. 1. As shown in the bottom panel, hi^^i^ 
values, as a function of ^3, for the two different reionization 
models, cannot be simply scaled by an overall normalization. 
Since is summed over the reionization visibility function, 
this is merely due to the difference in the reionization dura- 
tion. The same can be inferred through Fig. 2 where, as the 
redshift is increased, Ii curves shift to higher multipoles. 



lM^dMMj^b{M,z) 



(34) 



where M_ and M+ are the lower and upper limits of 
masses and 6(M, z) is the halo bias. We consider two 
possibihties for fluctuations in ^(r): (1) density inhomo- 
geneities and (2) patchiness. The bias factor related to 



the patchy model follows from Ref . [8] , while in the case 
of density inhomogeneities, we consider a constant bias 
factor of unity. This allows us to write bg{z) as 



bg{z) W l.O + bxAz) 



(35) 



where bx^ {z) = when Xe{z) = 1, such that we no longer 
consider patchiness of the reionization when the universe 
is fully reionized. The bias factor related to Xe fluctua- 
tions are described in Ref. [8] to which we refer the reader 
for further details; an interesting point related to this bias 
factor is that it is dominated by halos with temperature 
at the level of 10* K and above, where atomic cooling is 
expected and first objects form and subsequently reionize 
the universe. Since what enters in the cross-correlation is 
the bias factor times the growth of density perturbations, 
in this case, the product bxA^)G{^) is in fact a constant 
[25,26], as a function of redshift out to z of 20 or more 
given the rareness of halos at high redshifts. 

The source bias bs{z), for the three tracer field distri- 
butions in Fig. 1 again involves a similar description. In 
the case of our low redshift distribution — curve labeled 
(c) in Fig. 1 — we take bs{z) — 1, to be consistent with 
typical bias involved with galaxy fields. For the two high 
redshift ones, we calculate the bias following Eq. 34. For 
the broad high-z distribution, we take a bias factor con- 
sistent with the formation of first objects and set M_ 
to be the value corresponding to virial temperatures of 
lO^K and M_|_ 00. For the high-z narrow distribu- 
tion, we assume that neutral gas is only present in halos 
other than those that form first objects (and, thus, as- 
sumed to be reionized). Thus we take M_ and M+ 
the value corresponding to virial temperatures of lO^K; 
note that we have taken a simple description of source 
bias here for both halos containing first luminous objects 
and halos containing neutral Hydrogen. The situation 
is likely to be more complicated, but our main objective 
is to show that there is adequate signal-to-noise for a 
detection of the proposed polarization-temperature-high 
redshift cross-correlation. 

For models of extended reionization, the patchiness is 
more important. We illustrate this in Fig. 3 in terms of 
polarization power spectra. For reference, we show the 
primordial polarization anisotropy spectra generated at 
the last scattering, including the homogeneous rescatter- 
ing contribution at large angular scales. The secondary 
scattering contributions, due to inhomogeneities, can be 
written as [7,11] 



Ci = 



100 



"A 



99 



k= -^,r 
dA 



(36) 



where Q2^,(r) = J k'^dk/2Tr'^Q^°'){k,r). Note that the 
patchy polarization is generally higher [8] and may be 
the most important secondary polarization contribution 
related to scattering at arcminute scales. 

Fig. 3 also illustrates why the secondary polarization 
contributions may be undetectable from the power spec- 
trum alone; since the peak of the secondary contribution 



7 



is at the same angular scales as the primordial power 
spectrum peak, with an amplitude four-orders of mag- 
nitude smaller, the cosmic variance may not allow one 
to extract the secondary polarization contribution eas- 
ily. For comparison, in Fig. 3, we also show the expected 
noise power spectra related to two upcoming measure- 
ments of polarization involving Planck from space and a 
ground-based experiment, such as using a large- format 
polarization-sensitive bolometer array, with 1000 pixels, 
at the South Pole Telescope. We take one arcminute 
resolution, but limited to a sky-coverage of 4000 deg.^. 
Other possibilities include experiments such as QUAD 
[27], but with slightly lower resolution (at the level of 5 
arcminutes). 

Given that arcminute scale polarization information is 
correlated with the tracer field, one does not require all- 
sky maps of polarization or that of the tracer field. The 
required temperature information, however, comes from 
large angular scales, as we have shown in Fig. 2 with the 
integral of /*. The I value here corresponds to the tem- 
perature anisotropy multipole and what correlates with 
polarization fluctuations is the large angular scale tem- 
perature fluctuations. On the other hand, at the same 
time, arcminute-scale density fluctuations correlate with 
fluctuations in the visibility function and those respon- 
sible for secondary polarization signal. The large angu- 
lar scale temperature information, fortunately, is already 
available from the WMAP data to the limit allowed by 
the cosmic variance. While CMB polarization informa- 
tion, at arcminute scales, will soon be avilable both from 
ground and space, the extent to which the proposed study 
can be carried out will be limited from maps of the high 
redshift universe. 

As calculated in Eq. 25, the bispectrum is simply pro- 
portional to the quantity bi^i^ . This is an integral over the 
radial distance of the redshift distribution functions re- 
lated to the tracer fleld and electron scattering visibility, 
with their product weighted by the cross power spectrum 
of fluctuations between the two. To understand how this 
behaves, we plot fej^/a in Fig- 4 as a function of I2 for a 
given I'i (top panel) and as a function of ^3 for certain 
values of I2 (bottom panel) . These curves generally trace 
the integral shown in Fig. 2, but with an overall ampli- 
tude given by the strength of tracer field fiuctuations and 
captured with the index Z3. Here, wc show the difference 
between our two reionization models and using the broad 
high-z distribution for the tracer field. Since the reioniza- 
tion history related to model A is broader and expands 
to a higher redshift than model B, the function bi^i^ for 
model A, for any value of ^3, expands to higher values of 
h- As shown in Fig. 2, as the redshift is increased, I{l2) 
increases; since this redshift traces the reionization his- 
tory, we find that in histories where scattering happens 
at a higher redshift than a correspondingly lower one, the 
contribution to the bispectrum increases as a function of 

The same difference is partly responsible for the dif- 
ference in overall amplitude of bi^i^, but as a function 



Z3 with I2 fixed (lower panel of Fig. 4). Here, the over- 
all shape of the curve as a function of ^3, at a given I2, 
is determined by the clustering strength and mainly the 
shape of the cross-power spectrum between fluctuations 
in the scattering visibility function and the tracer field. 
Incidentally, we also note that the amplitude changes be- 
tween positive and negative values since the mode coup- 
ing integral, //, oscillates between positive and negative 
values. The transition from positive to negative values 
is simply, again, a refiection of the reionization history. 
Note that the difference in curves of Fig. 4 between the 
two reionization models cannot simply be considered as 
an overall scaling of the amplitude. The bispectra can, 
however, be considered as a scaling along the I2 axis for 
different reionization models. 



B. Signal-to-Noise Estimates 

In order to consider the extent to which the cross- 
correlations can be detected and studied to understand 
the reionization process, we estimate the signal-to-noise 
ratio for a detection of the bispectrum and smnmarize 
our results in Fig. 5. The signal-to-noise is calculated as 



E 

11I2I3 



^EE,tot ^CMB,tot ^SS ,tot 



12 



13 



(37) 



in the case of the E-mode related bispectrum while E 
B when the B-mode bispectrum is used. Here, (7*'*"* 
represents all contributions to the power spectrum of the 
ith field and we write 



^i,tot _ _|_ ^noise _|_ ^foreg 



(38) 



where (7"°""^ is the noise contribution and cf"^^^ is the 
confusing foreground contribution. We refer the reader 
to Ref. [12] for details on the derivation related to this 
signal-to-noise and the relation between bispectrum and 
various other statistics at the three-point level, such as 
the skewness. Instead of focusing on these statistics, 
which are all reduced forms of the bispectrum either in 
real-space or Fourier-space, we will focus here on the bis- 
pectrum directly and consider its detection. Note that 
the bispectrum, as discussed in Ref. [12], capture all in- 
formation at the three-point level and all other forms of 
statistics at this level capture only a reduced amount of 
information with a decrease in the signal-to-noise ratio 
depending on the exact statistic used and the shape of 
the bispectrum. While other statistics, such as skewness 
in real space is more easily measurable, we note that 
techniques now exist to construct the bispectrum reli- 
ably from CMB and large scale structure maps and have 
been successfully applied to understand the presence of 
non-Gaussianities in current data. Thus, we do not con- 
sider the measurement of the proposed bispectrum, in- 
volving CMB temperature, polarization and tracer-fields 
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FIG. 5. The cumulative signal-to-noise for the detection of the tcmporaturc-polarization-high redshift bispectrum, as a 
function of the multipole to which information is available for the high-redshift map of a tracer field of density fluctuations. In 
the left panel, we make use of model A to describe reionization history and plot signal-to-rioise ratios for Planck (all-sky; dotted 
lines), SPTpol (4000 sqr. degrees; dashed lines), and for a perfect experiment limited by cosmic variance only (all-sky; solid 
lines). Note that for different sky-coverages, signal-to- noise ratios scale as f'^-^. In each of these observational scenarios, we 
consider throe possibilities related to the tracer field labeled from the top to bottom as (a), (b) and (c), respectively following 
Fig. 1. The bottom dot-dashed line illustrates the decrease in the signal-to-noise ratio, compared to the to the top- most 
solid curve, with temperature information at multipoles greater than 25; the correlation is dominated by temperature at large 
angular scales, and this is already available from WMAP data. In the right-panel, we compare the extent to which reionization 
histories related to models A and B can be distinguished. Due to the lack of scattering at redshifts ^ 20 in model B, when 
compared to model A, we find that the cross- correlation is significantly reduced when using narrow z-band tracer fields. The 
planned upcoming studies at low radio frequencies involving the 21 cm background from neutral hydrogen prior to around the 
time of reionization may provide such narrow distributions of tracer fields. 



of higli-z universe, to be any more complicated than what 
is already achieved [30]. 

In Eq. 38, when describing noise, in the case of CMB 
temperature maps, we set cCMB^tot ^ ^cmb ^ ^noise^ 
where the noise contribution is related to WMAP and 
Planck data. Since for small values of l\ below few hun- 
dred, C^MB Qnoisc cTirrcnt data, we only need 
to consider cosmic variance. For polarization maps, we 
take C™'*°* = C™ + C;"°'^° + Cf , where C^^'^" is the 
polarization noise and Cf is the secondary contribution 
to polarization anisotropy power spectrum as shown in 
Fig. 3. In general, Cf^ » Cf , but C;™ is not signifi- 
cantly higher than C^°^^^ for upcoming experiments such 
as Planck such that one must account for noise properly 
beyond just the cosmic variance. For the tracer field, we 
primarily consider only the presence of cosmic variance 
such that cfS'*"* = C?^. This is due to the fact that 
we do not have detailed information related to detector 
and/or instrumental noise in certain high-z observations 
such as in the 21 cm background. 

We can roughly estimate noise in one case involving 
the planned measurement of IR fluctuations due to high- 
z luminous sources such as the population of first propto- 
galaxies. We include this noise when considering cross- 
correlation studies with the high-z broad distribution 
shown in Fig. 1, since such a broad redshift distribu- 



tion is expected for the high-z IR bright sources. In this 
case, we also include the shot-noise contribution to the 
tracer-field map associated with the finite density of these 
sources and based on models in Ref. [13]. Since the shot- 
noise generally peak at scales less than an arcminute, 
and the cross-correlation we are after is at angular scales 
of few arcminutes and more, we, however, find that this 
additional noise at small scales not to be important. 

As shown in Fig. 5 (left panel), if the reionization his- 
tory follows model A, cumulative signal-to-noise ratios 
range from about a few with Planck polarization maps 
to ^ lO'' using perfect maps of the arcminute polariza- 
tion field. These signal-to-noise ratios are substantial and 
are due to the fact that the correlation between tempera- 
ture and the temperature quadrupole, as well as the cor- 
relation between visibility fluctuations and maps of the 
high-z universe, is high. As shown in Fig. 2, if scattering 
were to be restricted to very low redshifts, z 1, then 
the correlation between temperature and the tempera- 
ture quadrupole is not significant at multipoles greater 
than a few. On the other hand, as one moves to redshifts 
~ 20, this correlation is broader and spans over few tens 
of multipoles. Considering this, the existence of the high 
signal-to-noise can be understood simply as following. If 
we rewrite Eq. 32 for the case with a narrow redshift dis- 
tribution, such as the case for high-z narrow correlation 
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FIG. 6. The cumulative signal-to-noise for the detection of 
the temperature-polarization-high redshift bispcctrum using 
maps of the B-mode polarization and assuming a perfect ex- 
periment and a tracer field given by the (a) distribution of 
Fig. 1. The top dot-dashed line is for the case where the lens- 
ing B-modo power is reduced with quadratic Icnsing statis- 
tics (see text for details). For comparison, we also show the 
maximum E-modo signal-to-noise ratio. For other distribu- 
tions considered, the signal-to-noise ratios scale in the same 
manner as shown in Fig. 5. In general, B-mode polarization 
maps provide a higher signal-to-noisc ratio than a map of the 
E-mode. Though the amplitude of the foreground reioniza- 
tion signal is same in both E- and B-mode maps, due to a 
significant decrease in the primordial signal which contribute 
to the overall variance, the cross-correlation study is better 
with B-modcs than E-modcs as the sample variance, relative 
to the foreground signal, is reduced. 

involving 21 cm fluctuations, say at a distance of r^, we 
have 

biuh,h^I^Ar,)Cf,'{r.), (39) 

where Cf^ {r^) is the cross-power spectrum between fluc- 
tuations in the visibility function and that of the tracer 

held, and is simply Cf (r,) - TgSsJcf^. Simi- 
larly, lf{rx) roughly scales as the cross-power spectrum 
between E-mode polarization and temperature, at r^, 

and we can write I^{rx) ^ tet^^ CfJ" . The signal- 
to-noise ratio scales as 

(I) ^Y.fi^^M-^^^^f§k^, (40) 

where /(Zi, Z2, Zs) ^ Q q ^ . The summation 

over Z3, l^Ci^ is simply the variance related to visibil- 
ity fluctuations and increases as the maximum value of I 
over which Z3 is summed. The ratio of Cj^^ / ((7™-*°*) be- 
haves on the signal-to-noise associated with polarization 



measurements, especially when li ^ I2. Thus, roughly, 
signal-to-noise in the correlation bispectrum, for each li 

and I2 mode, scales as /sVy ^l;T^gS'''gs(^/-'^)poi ^^nd is de- 
termined mostly by the resolution of the polarization 
map and the extent to which fluctuations in the visibility 
function is correlated with that of the tracer fleld. Since 
^ET ^ 1 ^^'^ '^gg is ^ 0-6 percent (when Z3 ^ lOOO)-'-, 
if r^g ^ 1, one can achieve, in principle, signal-to-noise 
values of order ~ few hundred with all-sky maps of po- 
larization and adequate resolution in polarization down 
to multipoles of 1000. This, of course, assumes no noise 
maps in both polarization and in the tracer field and 
just that the noise is limited by the cosmic variance. In 
Fig. 5, we include instrumental and observational noise 
which lead to a decrease in the signal-to-noisc ratio, for 
example, with Planck by several orders of magnitude. 

In Fig. 5 (right panel), we illustrate the difference be- 
tween expected signal-to-noise ratios for reionization his- 
tories involving models A and B. Here, we consider two 
tracer-fields involving the high-z broad distribution and 
the high-z narrow distribution. Fig. 1 shows the overlap 
between these distributions and the visibility function re- 
lated to the two reionization models. This overlap, or its 
non-existent, helps understand differences in the bispec- 
trum signal-to-noise ratios. In particular, we note the 
sharp reduction in the signal-to-noise ratio between the 
two models when using the narrow redshift distribution 
for tracer-field at high redshifts from model A to model 
B. This is due to the fact that most, if not all, rescatter- 
ing is concentrated at redshifts below 15 in the case of 
model B while the tracer fleld is a probe of structure at 
redshifts of order 20. 

The situation is also the same when one uses a tracer- 
field with a broad redshift-distribution, but limited to 
redshifts greater than 15. While this allows a mechanism 
to distinguish the extent to which two broadly different 
reionization models can be distinguished form one an- 
other, though they both give the same optical depth and 
essentially produce the same polarization signature, one 
can do more than this. Under the assumption that the 
reionization history follows an extended period, one can 
use tracer-fields with narrow band distributions to study 
the strength of the correlation as a function of redshift. 
This in return can be converted as a measurement of the 
visibility function as a function of redshift. While we have 
not considered such a possibility, due to the lack of in- 
formation related to how well narrow band distributions 
can be defined at redshifts of order 15, we emphasize that 
21 cm observations of the neutral Hydrogen content at 
high redshift may be the ideal way to approach such a 



'This can be understood l)ascd on the fact that the sec- 
ondary polarization due to visibility fluctuations shown in 
Fig. 3 is simply Ci ~ 3/lOOCf^QLs and since Qrms ~ 25 
/iK, is of order 4 x lO"'' or cigg ~ 0.6% 
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reconstruction; Due to the line emission, redshift ranges 
can be a priori defined at observational wavelengths and 
such observations, when combined with CMB, may pro- 
vide an interesting possibility for detailed studies related 
to reionization. 

Instead of E-modes, as mentioned, we can also consider 
the B-mode for a cross-correlation study (Fig. 6). While 
there is no large angular scale signal related to rescatter- 
ing in the B-mode map, associated with density fluctu- 
ations, at tens of arcminute scales and below, a B-mode 
signal exists due to modulations of the visibility function 
by density inhomogeneities in the ionized electron distri- 
bution. For the B-mode bispectrum, to describe noise, 
we take cfj^''^"* = Cll^ + Cf^''^'' + Cf and set Cj^^ as that 
due to gravitational lensing conversion power from E- to 
B-modes, as shown in Fig. 3, and ignore the presence of 
primordial gravity waves or tensor modes. This is a safe 
assumption given that such primordial contributions are 
expected to peak at angular scales of few degrees, while 
the proposed study involve tens of arcminute scale polar- 
ization. 

Since the ratio of Cf^'"^/(C;^^) is higher than that re- 
lated to E-modes, one generally obtains a higher signal- 
to-noise ratio with an arcminute scale B-mode map when 
compared to a map of E-modes. The same is expected 
since the primordial fluctuations are smaller in B-modes 
than in E-modes, while the secondary contribution re- 
mains the same, the confusion generated by primary fluc- 
tuations for the cross-correlation is reduced in the case 
of B-modes relative to the same in E-modes. This can 
also be stated in terms of the reduced sample variance. 
There is another advantage with the use of B-mode in- 
formation. Since the main confusion at arcminute scales 
is related to that associated with gravitational lensing 
confusion, one can improve the signal-to-noise ratios for 
the bispectrum measurement significantly, when the B- 
mode map is a priori cleaned with quadratic statistics 
or improved likelihood techniques [29] . These techniques 
lower the B-mode power related to lensing by more than 
an order of magnitude and reduces the sample variance 
contribution to the noise associated with the bispectrum 
by a similar number. 

C. Contaminant Correlations 

Using Figs. 5 and 6, we have argued that there is suffi- 
cient signal-to-noise for the observational measurement of 
the proposed correlation in the form of a bispectrum. It is 
also useful to consider if this bispectrum can be confused 
with other sources of non-Gaussianities in the CMB tem- 
perature and polarization anisotropics when correlated 
in the same manner with a map of the high redshift uni- 
verse. The only possibility is related to the weak lensing 
effect on CMB anisotropies [28]. Essentially, weak lens- 
ing deflections of CMB photons lead to a correction to 
CMB anisotropy that depends on the temperature gra- 




FIG. 7. The cumulative signal-to-noise for the detection 

of the tcmporaturc-polarization-high redshift bispectrum, re- 
lated to the gravitational lensing effect on CMB tempera- 
ture and polarization anisotropies. The uon-Gaussiau sig- 
nal is generated due to the fact that the high-redshift tracer 
field may correlate with the lensing deflection potential that is 
modifying the spatial distribution of CMB fluctuations. Wo 
show the signal-to-noise ratio for the tracer held involving 
the high-z broad distribution and for a perfect, Planck and 
SPTpol like ground-based experiment. The top curves are for 
the case where all information from the temperature map is 
included while the bottom curve, for a perfect polarization 
measurement, is for the case where we restrict temperature 
information to a multipole of 100 only. While the reionizar 
tion related bispectrum is not affected by such a cutoff, due to 
the fact that relevant scales related to the correlation comes 
from large angular scale temperature anisotropy, in the case 
of lensing, this results in a significant reduction of this sig- 
nal-to-noise. The dot-dashed line shows the signal-to-noise 
ratio associated with a bispectrum measurement related to 
lensing using the B-modo map, in addition to temperature 
and traccr-flcld information. Note the difference in scale of 
the ordinate between Figs. 5 and 6 and this plot. 

dient and the deflection angle. While the temperature 
gradient correlates with the E-mode polarization map, 
the deflection angle can correlate with a map of fluctua- 
tions in the high-z universe, and generate a non-Gaussian 
signal related to the bispectrum. 

We write the bispectrum related to this cross- 
correlation as 

^ 1 / (2/l + l)(2/2 + l)(2^3 + l) 

hhh 2 V 47r 

{(20 J^)n^l,^2,/3)C«^Cf 

+ (0 '^)nh,h,hK'ctf}, 

(41) 
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where F{hMM) = Mh + \) + h{h + I) - hih + \)] 
Cf^ is the cross-power spectrum between lensing poten- 
tials and the tracer field and is given by 



where 



t:./_ ^ dA{ro - r) 



W^{k,r) = -m^[^] F{r)- 

k J dA{r)dA{ro) 



(42) 



(43) 



Similarly, lensing induces a bispcctrum related to B- 
modes combined with temperature information and a 
tracer-field. We can write this bispectrum as 



BOS 
I1I2I3 



i /(2;i + i)(2;2 + i)(2;3 + i) 



(44) 



In general, for high-z distributions of the tracer-field, 
W^{r) and do not overlap significantly. This results 
in a reduction in the cross-correlation between the two. 
Moreover, since the reionized polarization anisotropy re- 
lated bispectrum can be constructed with Qi < 100 or 
so without a loss of signal-to-noise, the importance of 
the lensing bispectrum can be significantly reduced since 
limiting the lensing bispectrum to I2 values less than 100, 
leads to a substantial reduction in the signal-to-noise as- 
sociated with it. Even with a perfect experiment, and I2 
multipole taking same values as Z3 out to 1000 or more, 
the signal-to-noise ratios arc not greater than few tens 
[28] (see. Fig. 7). With I2 reduced to a smaller value be- 
low 100, this signal-to-noise ratio decreases substantially 
below a few, at most. Given that the polarization fluc- 
tuations related to bispectrum has substantially higher 
signal-to-noise ratios, it is unlikely that the lensing gen- 
erated non-Gaussianities are a major source of concern 
for the proposed study related to reionization. 



While the fluctuation amplitude, relative to dominant 
primordial fluctuations, is small, we suggest that a cross- 
correlation between arcminute scale CMB polarization 
data and a tracer field of the high redshift universe, 
such as through fluctuations captured by the 21 cm neu- 
tral Hydrogen background or those in the infrared back- 
ground related to first proto-galaxies, may allow one to 
study additional details related to reionization. This in- 
cludes a possibility to recover the visibility function re- 
lated to electron scattering. For this purpose, we discuss 
an optimized higher-order correlation measurement, in 
the form of a three-point function, involving information 
from large angular scale CMB temperature anisotropies 
in addition to arcminute scale polarization data. The 
proposed bispectrum can be measured with a substan- 
tial signal-to-noise ratio with lens-cleaned B-mode maps 
where the confusion related to primordial anisotropies 
and the arcminute scale reionization-related polarization 
is the least; this can be compared to other sources of non- 
Gaussianity in CMB data, such as due to gravitational 
lensing, which with the same combination of tempera- 
ture, polarization, and the tracer field is at the level of 
at most unity in terms of the signal-to-noise ratio. 

Given that the arcminute scale polarization fluctua- 
tions are correlated with fluctuations at high redshift at 
same angular scales, for a reliable measurement of the 
three point correlation function suggested here, one does 
not require all-sky maps of CMB polarization or that 
of the tracer field. This is helpful since one can ob- 
tain high signal-to-noise maps of the polarization from 
ground-based experiments, and in certain cases fluctua- 
tion data on certain tracers, by concentrating on limited 
sky-area instead of the whole sky. The required infor- 
mation from temperature anisotropies is at large angular 
scales and this information is already available to the 
limit allowed by the cosmic variance with WMAP data 
and is expected to improve with Planck in terms of the 
foreground confusion and separation. A study such as 
the one proposed may allow one to establish the epoch 
when CMB polarization related to reionization is gener- 
ated and to address if the universe was reionized once or 
twice. 



IV. SUMMARY 



ACKNOWLEDGMENTS 



The free-electron population during the reionized 
epoch rescatters CMB temperature quadrupole and gen- 
erates a now well-known polarization signal at large an- 
gular scales. While this contribution has been detected in 
the temperature-polarization cross power spectrum from 
WMAP data, due to the large cosmic variance associ- 
ated with anisotropy measurements at relevant scales, 
only limited information related to reionization, such as 
the optical depth to electron scattering, can be extracted. 
The inhomogeneities in the free-electron population lead 
to an additional secondary polarization anisotropy con- 
tribution at arcminute scales. 



This work was supported in part by DoE DE-FG03- 

92-ER40701 and a senior research fellowship from the 
Sherman Fairchild foundation. 



[1] M. Zaldarriaga, Phys. Rev. D 55, 1822 (1997) 
[arXiv:astro-ph/9608050] . 



12 



[2] C. L. Bennett et al, Astrophys. J. Suppl. 148, 1 (2003) 

[arXiv:astro-ph/0302207]. 
[3] A. Kogut et al, Astrophys. J. Suppl. 148, 161 (2003) 

[arXiv:astro-ph/0302213]. 
[4] R. Cen, Astrophys. J. 591, 12 (2003) [arXiviastro- 

ph/0210473]. 

[5] M. Kaphnghat, M. Chu, Z. Haiman, G. Holder, L. Knox 
and C. Skordis, Astrophys. J. 583, 24 (2003); G. Holder, 
Z. Haiman, M. Kaplinghat and L. Knox, Astrophys. J. 
595, 13 (2003); W. Hu and G. P. Holder, Phys. Rev. D 
68, 023001 (2003). 

[6] e.g., see the review by R. Barkana and A. Loeb, Phys. 
Rep. 340, 291 (2001). 

[7] W. Hu, Astrophys. J. 529, 12 (2000) [arXiviastro- 
ph/9907103]. 

[8] M. G. Santos, A. Cooray, Z. Haiman, L. Knox and 

C. P. Ma, arXiv:astro-ph/0305471. 
[9] S. P. Oh, arXiv:astro-ph/0005262. 
[10] S. H. Hansen and Z. Haiman, arXiv:astro-ph/0305126; 

X. Chen and M. Kamionkowski, arXiv:astro-ph/0310473; 

S. Kasuya, M. Kawasaki and N. Sugiyama, arXivrastro- 

ph/0309434. 

[11] D. Baumann, A. Cooray and M. Kamionkowski, New As- 

tron. 8, 565 (2003) [arXiv:astro-ph/0208511]. 
[12] e.g., in the case of CMB temperature alone, A. Cooray 

and W. Hu, Astrophys. J. 534, 533 (2000). 
[13] A. Cooray, J. J. Bock, B. Keating, A. E. Lange and 

T. Matsumoto, arXiv:astro-ph/0308407. 
[14] P. Tozzi, P. Madau, A. Meiksin and M. J. Rees, 

arXiv:astro-ph/9903139. 
[15] X. Fan, V. K. Narayanan, M. A. Strauss, R. L. White, 

R. H. Becker, L. Pentericci and H. W. Rix, arXivrastro- 

ph/0111184; R. L. White, R. H. Becker, X. H. Fan and 

M. A. Strauss, arXiv:astro-ph/0303476. 
[16] R. Cen, Astrophys. J. 591, L5 (2003); Z. Haiman and 

G. P. Holder, Astrophys. J. 595, 1 (2003); J. S. B. Wyithe 

and A. Loeb, arXiv:astro-ph/0302297. 
[17] W. H. Press and P. Schechter, Astrophys. J. 187, 425 

(1974). 

[18] X. Chen, A. Cooray, N. Yoshida and N. Sugiyama, 

arXiv:astro-ph/0309116. 
[19] ,J. N. Goldberg, J. Math. Phys. 7, 863 (1967) 
[20] M. Zaldarriaga and U. Scljak, Phys. Rev. D 55, 1830 

(1997); M. Kamionkowski, A. Kosowsky and A. Stebbins, 

Phys. Rev. D 55, 7368 (1997) [arXiv:astro-ph/9611125]. 
[21] R. K. Sachs and A. M. Wolfe, Astrophys. J. 147, 73 

(1967). 

[22] A. Cooray and R. Shcth, Phys. Rept. 372, 1 (2002) 
[arXiv:astro-ph/0206508]. 

[23] D. J. Eisenstein and W. Hu, Astrophys. J. 511, 5 (1997) 
[arXiv:astro-ph/9710252]. 

[24] E. F. Bunn and M. J. White, Astrophys. J. 480, 6 (1997) 
[arXiv:astro-ph/9607060]. 

[25] P. J. E. Peebles, The Large Scale Structure of the Uni- 
verse (Princeton: Princeton Univ. Press) (1980). 

[26] S. P. Oh, A. Cooray and M. Kamionkowski, Mon. 
Not. Roy. Astron. Soc. 342, L20 (2003) [arXiviastro- 
ph/0303007]. 

[27] M. Bowden et al, arXiv:astro-ph/0309610. 

[28] W. Hu, Phys. Rev. D 62, 043007 (2000) [arXiv:astro- 



ph/0001303]. 

[29] T. Okamoto and W. Hu, Phys. Rev. D 67, 083002 
(2003) [arXiv:astro-ph/030103li; M. Kcsden, A. Cooray 
and M. Kamionkowski, Phys. Rev. D 67, 123507 (2003) 
[arXiv:astro-ph/0302536]; C. M. Hirata and U. Sel- 
jak, Phys. Rev. D 68, 083002 (2003) [arXiv:astro- 
ph/0306354]. 

[30] E. Komatsu, B. D. Wandeh, D. N. Spergel, A. J. Ban- 
day and K. M. Gorski, Astrophys. J. 566, 19 (2002) 
[arXiv:astro-ph/0107605]; M. G. Santos et al, Phys. Rev. 
Lett. 88, 241302 (2002) [arXiv:astro-ph/0107588]. 



13 



